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Abstract 

Following our previous work wherein the leading order effective action was computed in 
the covariant effective field theory of gravity, here we specialize the effective action to the 
FRW spacetime and obtain the effective Friedmann equations. In particular, we focus our at¬ 
tention on studying the cosmological implications of the non-local terms when each of them is 
combined with the Einstein-Hilbert action. We obtain both analytical and iterative solutions 
to the effective background equations in all the cases and also briefly comment on the consis¬ 
tency between the iterative and numerical solutions whenever possible. We find that among 
all the non-local terms, the imprints induced by are very significant. Interpreting 

these corrections as an effective dark energy component characterized by an equation of state 
parameter, we find that the R^R correction can indeed lead to an accelerated expansion of 
the universe at the present epoch even in the absence of a cosmological constant. We briefly 
discuss some phenomenological consequences of our results. 


1 Introduction 

Over the last few decades, cosmology has become a very precise science and it is indeed the right 
moment to move from phenomenology to theory, since now we finally have the possibility to falsify 
the proposed modifications to Einstein’s theory of general relativity (GR). While a large class of 
such modifications suggested recently are purely phenomenological, it is possible to predict these 
corrections from first principles within the framework of effective field theory (EFT), even in the 
absence of a complete theory of quantum gravity EEIEIII!. When applied to gravity, the EFT 
framework not only allows us to consider it together with the other forces of nature but it also 
manages to make distinct predictions. When the EFT of gravity is constructed in a complete 
covariant approach, as we have done in our previous paper [5], it is indeed possible to obtain 
these modifications to GR on a general background and in particular, in a cosmological setting. 

One of the major problems in cosmology is to explain the accelerated expansion of our universe 
at the present epoch, which is confirmed by different cosmological probes and has motivated a 
lot of activity in order to provide a satisfactory theoretical explanation for this observational fact 
[HIE]. Among the plethora of possibilities, the cosmological constant A is unarguably the simplest 
choice, however, the typical problems associated with it e.g. its enormously small value and its 
dominance only today has led to the exploration of other potentially viable models [giaiini[ii]. 
It is well known that cosmic acceleration within the framework of GR can be achieved either 
by adding an ‘exotic’ matter component to GR or by the modification of GR itself. As a result, 
these alternative scenarios can be broadly classified into two main categories: i) scalar field models 
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nails] and ii) modified gravity models |14 I I15 | [T6]. In the former case, either a minimally coupled 
scalar field with a suitably flat potential is added to Einstein’s gravity or it is more entwined with 
gravity as, for instance, in the scalar-tensor theories m 118) . In the latter, the modifications 
to gravity could be local as in f{R) models as well as non-local type. Both of these broad 
possibilities have been studied to a great extent in the literature, for some recent reviews, see 

naEoi. 

As outlined earlier, the EFT of gravity leads to consistent (local and non-local) modifications 
of GR which could turn out to very relevant in different contexts and some of them can also 
explain dark energy. The distinction we want to make is that EFT induced corrections are 
exactly determined by computation, leaving no space for further tuning of the parameters. By 
adapting to a Friedman-Robertson-Walker (FRW) background and using the leading order (LO) 
effective action computed in [5], we start the investigation of the cosmological implications of 
the covariant EFT of gravity. In particular, we analyze the influence of the leading low energy 
non-local corrections to the background expansion and understand their role in the evolution of 
the universe at late times. 

The imprints of non-local terms have been considered in cosmology in many different contexts, 
for instance, as candidates of dark energy to explain the present acceleration of the universe 
[iiiiisiisaisaissiiMiisziiMiisaisniiaiEg, as a mechanism for screening of the cosmological 
constant |33[ 134) . to obtain accelerating solutions by means of a time delay in the Friedmann 
equations [35], for a resolution of the cosmological singularity IMIET], in the Newtonian limit 
|38) [39] and also to obtain the spherically symmetric solutions @01 HD [12]. Accelerating solutions 
have also been obtained in a general class of non-local higher derivative gravity [431 Hll Il5] . While 
this indicates that the non-local terms indeed lead to very rich and interesting phenomenology, 
they have mostly been studied as modifications to GR, and also argued as being consistent in 
some cases but have not been obtained in an effective covariant framework. 

This paper is organized as follows: in the following section, we quickly recap the EFT of 
gravity, write down the effective action to the quadratic order in curvature in the Weyl basis and 
specialize it to the FRW background. In section we study the classical theory as well as the 
leading local correction in the form of B? and obtain the analytical and numerical solutions. In 
section we study the leading non-local terms that arise at the second order in curvature and 
obtain the analytical, iterative and numerical solutions for all the cases. In sectionwe compare 
the these corrections in terms of the equation of state parameter for dark energy. Finally, in 
section we summarize our results and conclude with some discussions and a future outlook. 
In the appendix, we discuss some subtleties associated with the Green’s function in the FRW 
spacetime and outline how to obtain the kernel for the logarithmic distribution. 

2 Effective action from the covariant EFT of gravity 

In the previous paper [5], we had computed the leading order effective action to the second 
order in the curvature expansion within the covariant EFT of gravity. The final result for the 
gravitational part of the effective action is given by 

-J d‘x^RT(^'^R- + ( 1 ) 

where G, A, ^ and A are phenomenological parameters which must be fixed by the experiments or 
observations while the structure functions R and G are completely determined once the matter 
content of the theory is specified. The general form of these structure functions can be found in 
|5]. In order to understand the implications of this effective action in a cosmological setting, we 
work in a (3 + 1)-dimensional, spatially flat, FRW universe described by the line element 

= —dt^ + a^{t)dy? , ( 2 ) 
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Table 1: Coefficients of the non-local terms for various spins taken from [5]. Spin 0 fields are 
conformally coupled for x = 1 and minimally coupled for x = 0. Spins ^ and 1 are minimally 
coupled while the last column for spin 2, shaded in grey, are indeed the predictions of the EFT 
of gravity. All the numbers in the table must be divided by (dvr)^. In case, one considers N 
particles of a given matter species, the numbers in the first three columns should be multiplied 
by Nq,Ni /2 3 -nd A^i, respectively while the last column remains unchanged. 


where t is the cosmic time and a{t) is the scale factor. For this metric the Weyl tensor vanishes 
identically and therefore, the Weyl part of the effective action 0 will not contribute to the 
background equations of motion (EOM). However, this is not true when dealing with cosmological 
perturbations and the contributions due to these terms must be taken into account which could 
lead to very distinct signatures in the cosmological observables. As discussed in [5], at third order 
in the curvature expansion, the contributions from the conformal anomaly will also be present 
which can have interesting cosmological implications (see for instance [Ml HZ] and references 
therein) but we will not discuss them here. 

In this paper, as we are only interested in studying the background evolution, we can safely 
neglect the Weyl terms in ([^ and the total effective action including the matter action then 
simplifies to 

r = (i? - 2 A) - ^ y j + (3) 


where the quadratic R^ term is a local term while the structure function R is non-local in the 
low energy massless limit, i.e. m? <C —and has the following form 


R 



, — □ ^ — □ 

alog^ +/3— + 7—log^ + (5 

—U —U 



(4) 


The coefficients a, (3 ,7 and 5 are indeed predictions of the EFT of gravity which ultimately 
depend only on the field content of the theory and are listed in Table In the case of scalars and 
fermions, m is the effective mass of the particle. In the case of photons only the massless limit can 
be considered w? —)• i.e. only the first logarithmic term survives, while in the case of gravity, 
the effective mass is = —2A + 1Q'kGV{v). In the absence of a cosmological constant and/or 
an effective potential, we add an IR regularization mass p? as in the case of massless fields. 

In a cosmological context, the first logarithmic term, the only one present in the massless 
limit, has been studied in |21l 148] and in the framework of EFT in |49j . The successive ^ term 
is part of the class of non-local modifications of gravity first proposed in [ 22 ] while the last term 
in Q has been proposed and studied in |28l I31j . Apart from the leading logarithmic term, the 
other non-local terms were considered at a phenomenological level argued as consistent non-local 
modifications of gravity, but as shown in [5], all these terms arise in the framework of the EFT 
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of gravity and all of them must be taken into account together to construct a viable scenario of 
the universe. As mentioned earlier, the coefficients of these terms are not free parameters and 
one must tune the matter content of the theory for a given term to be the leading one. Below, 
we discuss the conditions under which this happens. 

• For the case of massless matter (without the graviton), only the scalar leading logarithmic 
term is present as shown in the Table However, it vanishes in the conformal case x = 1. 
Therefore, there is no correction from the structure function T and the leading term in the 
action will be the one induced by the conformal anomaly. 

• The case of massive matter (without the graviton) with x = 1 again eliminates the logarith¬ 
mic term thereby making the leading contribution being Also, if one takes A^i /2 = 2 A^i, 
all the logarithmic contributions vanish i.e. a = 7 = 0 . 

• One can find a condition for a theory (without the graviton) such that a particular coefficient 
vanishes. For instance, for x = 0 12Ao = 3 A^i /2 = 4A^i, one gets /3 = 0 while for 
A^o = 2 A^i /2 = 2A^i, one gets 7 = 0 . However, it is not possible to have both (3 and 7 vanish 
simultaneously. Also, for x = 1; one can obtain /3 = 0 for 15A"o = 6 A^i /2 = lOiVi and 7 = 0 
for Nii 2 = 2A^i. Moreover, for x = 1; it is instead possible to kill both (3 and 7 by choosing 
Nq = 3 A'i /2 = 6 A^i such that the leading correction is given by 

• In the presence of the graviton, it is not possible to make the leading logarithmic corrections 
go away, i.e. a = 0 , as the contributions due to both the scalars and gravitons are positive. 

• Including the graviton, it is possible to get (3 = 0 hy choosing 165A^o = llA^i /2 = ISA^i for 

X = 0 and by choosing SSNq = llA^i /2 = ISA^i for x = 1- One can also attain 7 = 0 by 
choosing 25Ao = A^i /2 = SOA"! for x = 0 2Ni/2 = 53Ai for x = 1- Furthermore, it is 

possible to make both f3 and 7 vanish simultaneously only for x = 1- 

Along the same logic presented above, one can also find a condition such that the correction 
term with 6 vanishes with or without the graviton. It is also evident from this analysis that when 
the graviton contribution is taken into account, it is not possible to make either 7 ^ or ^ the 
leading correction since the logarithmic term will always be present. 

In this paper, even if we have just shown that all the non-local terms should be considered 
together, we start by studying the solutions corresponding to different cases where each of these 
non-local terms is combined with the Einstein-Hilbert term separately and later, compare them 
with each other. For each case, we write down the covariant form of the field equations and 
obtain their solutions. In general, when a correction term is present in a theory along with the 
classical counterpart, the modified Einstein’s equation can be written as 

G^lu + = SirGTfj^i,, (5) 

where AG^iy is the correction term arising from the local or the non-local terms. Note that, the 
correction term AG^y is covariantly conserved. 

3 Local theory 

In this section, we first review the background solutions in the classical theory described by 
the Einstein-Hilbert action (with the cosmological constant A) and matter. Later we study the 
solutions for E? gravity which is the leading local correction to the classical theory within the 
EFT of gravity. 
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3.1 Classical theory 

At the classical level, the theory is described by the Einstein-Hilbert action (with the cosmological 
constant A) 

( 6 ) 

(7) 

( 8 ) 


r = ^ J d^x^g [R - 2A) + 

The Friedmann EOM can be obtained by varying this action 

2 ^ SttG 

= - H- p 

3 3 ^ 


a 
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A AttG 


(P + 3p), 


where H = d/a is the Hubble parameter and p and p indicate the energy density and pressure 
of the matter content, respectively. Conservation of the energy-momentum tensor leads to the 
continuity equation 

p + 3H{l + w)p = ^, (9) 

where w = p/pis the equation of state parameter of the matter content. Note that only two among 
the equations 0, 0 and 0 are independent. In a universe dominated by the cosmological 
constant with w = —1, the solutions to these equations lead to a de Sitter expansion, for which 
the scale factor is given by 


aA(t) = eV 3 . 

(10) 

1/3) with pr oc a~^, we get 


Mt)=(J , 

(11) 


while for the case of matter or dust {w = 0) with pm oc a one finds 

2/3 


— { 


(a' 


( 12 ) 


where the normalization is fixed by imposing a(to) = 1- In the later sections we will evaluate the 
leading order corrections on these background solutions. 


3.2 gravity 

As explained earlier, the local quadratic term appears at the leading order in the EFT of 
gravity 0. This is one of the most interesting and well studied example of a local higher order 
correction to Einstein’s theory which was originally studied in m- The action in this case is 


r = 


1 


lOvrG 


d^xy/^ {R-2K)--j d^x^R^ + 5, 


(13) 


where is a free parameter which must be fixed by, for instance, the condition that inflation is 
driven by this term. However, one can also think of the usual situation where inflation is induced 
by the matter sector and the R? term is not present at all. Although it is possible to study 


the dynamics of the action (13) in the Einstein frame by doing a conformal transformation, we 


rather choose to work in the Jordan frame as it will be useful for the case when this local term 
is combined with the non-local terms to study the dynamics of the full effective action 0. 

To obtain the correction term to the Einstein’s equations, we make a variation of the second 


term in the action (13) using the following relations 




6R = - V - > 


(14) 
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where = 5g^y. The resulting expression becomes 




= 2 


{R^lu - R-2 (V^Vj. - g^,yU) R. 


By inserting the FRW metric in this equation and using the relations 

Roo = -3(i^ + H‘^) R = 6{H + 2H^), 

in ^ leads to the modified Friedmann EOM for R^ gravity 


(2HH + 6H^H - hA = - + —. 
C \ / 3 3 ^ 


(15) 


(16) 


(17) 


which can be immediately compared with |50L EH- Note that, the inclusion of the correction 
term makes this differential equation of third order for the scale factor instead of first order in its 
absence. 

In the absence of matter, the EOM for R^ gravity is given by 


aH'^ + 2HH + QH'^H - H‘^ = 0 


a = — 




96ttG 


(18) 


This is a second order non linear differential equation which is not exactly solvable. In the 
following, we will first hnd exact analytical solutions for the pure R^ gravity and later obtain 
approximate and numerical solutions to this equation. It is useful to note that the classical theory 
is recovered as a —)• oo. 

In the case of pure R^ gravity (a = 0) we can find an analytical solution which for particular 
choices of the integration constants reduces to both de Sitter and radiation. We want to solve 


2HH + 6H‘^H -7/2 = 0. 


(19) 


Clearly H = H{tQ) is the constant solution to this equation which corresponds to the de Sitter 
space thus pure R^ gravity inflates, as is well known. To find the other solutions, we first shift 
to a new variable Z = H/ AH to find 


dZ = -SVHdH, 

which can be easily integrated to give the following implicit solution 

^ f dH 

t + C2 = 


CiVH - 2H^ 


( 20 ) 


( 21 ) 


depending on the two integration constants Ci and C 2 - The integral in (21) can be performed 


exactly but the resulting expression is not very illuminating apart from the case Ci = 6*2 = 0 , 
yielding 


H = 


1 

Yt 


( 22 ) 


which is indeed the radiation solution (Q. The fact that radiation is an exact solution of pure 
i?2 gravity is expected since R = 0 for radiation in FRW and any action with more than one 
power of R will lead to EOM that admits such a solution. The de Sitter solution is recovered in 
the limit C 2 —>• —00 when Ci 7 ^ 0 . 

As mentioned earlier, we can not solve exactly R? gravity in the absence of matter but we can 


try to understand the solutions in different regimes when certain terms in (18) can be neglected 


as compared to the others. If at early times H is nearly constant, we can ignore the higher 
order time derivative terms 2HH and as compared to 6H^H in (18), one finds the quasi 
de Sitter regime with the solution given by H{t) ~ H{to) — — to). On the contrary, when 

H is time dependent, the cubic term 6H^H can be dropped, leading to an oscillatory solution 
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Figure 1: On the left, the time evolution of the Hubble parameter H is plotted for gravity 
without matter. The dashed black line is the full numerical solution to (18). The dashed blue 
line indicates the constant solution for H normalized to H(to) while the solid blue line indicates 
the quasi de Sitter solution described in the text. The solid orange line represents the damped 
oscillating solution at late times and the dashed orange line indicates the slope ^ as in the case 
of matter. On the right, the time evolution of H for gravity including A is plotted as shown 
by the solid blue line. The presence of A, as indicated by the dashed blue line, suppresses the 
damped oscillating solution, leading to a late time de Sitter phase. 


. As shown in Figure these oscillations are indeed present when equation (18) is integrated 
numerically. In the presence of matter, it is easy to realize that R? gravity with cosmological 
constant or radiation is exactly solved by de Sitter and radiation solutions, respectively. 

We now turn our attention to solve R? gravity numerically. We first solve it without any 
matter content and then solve it in the presence of A as well as with radiation and matter. For 
these cases, we solve the EOM ( pT[ ) together with the continuity equation ([^. For R? gravity 
without matter, we find that the exact numerical solution shows a quasi de Sitter phase followed by 
a decaying oscillatory behavior and is well corroborated by the approximate solutions in different 
regimes discussed earlier, as shown in Figure The effect of adding the cosmological constant 
is to damp the oscillations at late times before reentering the de Sitter phase. In Figure we 
show the solutions for R? gravity with radiation as well as with matter and radiation. In both 
cases, the oscillatory phase damps out and the solution enters into a radiation phase. Note that 
the global behavior of these solutions is qualitatively very interesting from a phenomenological 
point of view. 


4 Non-local theory 

In this section, we shall combine the Einstein-Hilbert action with different non-local terms present 
in study the solutions of the relative effective Friedmann EOM and discuss their relevance 
for the late time evolution of the universe. 

4.1 i? log gravity 

As discussed in our previous paper [5], the leading non-local correction appearing at the second 
order in curvatures in the EFT of gravity is given by the logarithmic term which has been studied 
in [211 EH] and in the framework of EFT in [49j . In this section, we will study the effects of the 
non-local action of the form 

Inl = [ Rlog-^R. (23) 
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Figure 2: On the left, the time evolution of the Hubble parameter H is plotted for gravity 
together with radiation. After the oscillatory solution damps out, both H and the radiation 
energy density pr asymptote to ^ and as shown by the dashed blue and dashed orange lines, 
respectively. In the right panel, the time evolution of H is plotted for R? gravity with both 
radiation and matter present. Also for this case, as expected, after the end of the oscillatory 
phase, H, pr and the matter energy density pm nicely settle to their asymptotic behaviors, as 
indicated by the various dashed lines, respectively. 


As discussed in section this correction term is only generated by minimally coupled scalars 
and gravitons while it vanishes for fermions, photons and conformally coupled scalars. Before 
proceeding, we first clarify the precise meaning of the logarithmic term in (23) which is given by 
the following integral representation 


iog^= r 
Ai Jo 


ds 


1 


1 


p,^ + s —□ 


(24) 


This dehnition allows us to express the logarithm in terms of the Green’s function of the operator 
—D+s, from which it will inherit the boundary conditions. We now provide a covariant derivation 
of the EOM for this case. Upon variation, the action in (23) leads to 
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SInl = / < 6^/^Rlog -^R + 2^/^^Rlog -^R + ^/^R6log -^R 
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term, we vary (24) to arrive at 


The first two terms can be easily dealt with using the variations given in (14) while for the last 

(26) 
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where the variation of □ acting on a scalar is given by 

50(1) = - . 

Using these relations, it is now easy to extract the modification to the Einstein’s tensor 
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When dealing with non-local corrections of this type, it turns out to be very useful to introduce 
auxiliary fields to localize them at the expense of increasing the number of equations. In this 
particular case, we define L = log and Us = zi^^R to obtain 


IGvrGa 


L 

(-□ + s)Us 


2 { Rij.u - -^QiiuR ] L-2 L 


ds 


—+ -g/iu'^'^UsVaUs + -gf^pUgOUs 


poo 

r 1 1 

/ ds 

, R Us 

Jo 



= R. 


Note that, while the first line of (28) is equivalent to (15) after the substitution L 


(28) 

(29) 

(30) 
R, the 


integral term in the second line appears due to the variation of the logarithm. 


We are interested in obtaining the correction term (28) in the FRW background. The time¬ 
time component of the correction term is 


1 /1 \ 1 


-ut + UsUs + SHUsUs 


(31) 


where the first two terms are evaluated using the following expression for the field L 

L{t)= [ dt'a^{t')Ll^"^{t-t')R{t'), (32) 


and the kernel is the retarded representation of the operator log in the FRW background. 

Ignoring the time dependence of the scale factor it can be written as 


such that the FRW kernel is just the flat space kernel 

9{t — t' — e) 


LyHt - t') = -2 lim 

e-s>0 


t - t' 


+ 5{t - t') log/re 


(33) 


(34) 


except the time dependent normalization. The details of this derivation can be found in appendix 
[a} For the moment, we also neglect the last integral term in equation (31), which appears due 
to the variation of the logarithmic term. Note that, this term is a higher order term (in time 
derivatives) as compared to the first two terms. The resulting EOM by including this term will be 
a non-trivial integro-differential equation. For our purpose, to access these quantum corrections 
and to keep the calculation tractable as much as possible, we have used the approximation 
wherein the higher order time derivative terms can be ignored. We now evaluate the first two 


terms by inserting (32) together with (33) and (34) into (31) which leads to the following effective 
Friedmann EOM 

+ IQ-kGu (^R{t) + a~^{t) J dt'a^{t')R{t') Ly^{t — t') 

— 327rGaH{t) a~2[t) j dt' Ly^{t — t') ^ {t')R{t')j = ^ ^^ p(t), (35) 


which can be immediately compared with [39] . Although the correction term appears non-local 
in time due to the retarded boundary conditions imposed, it is unambiguous and perfectly causal. 


However, the effective EOM (35) is an integro-differential equation and as such it is generally 
difficult to solve. 
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Radiation is an exact but trivial solution since for R = 0 the correction term in (35) vanishes 
(as in the case of gravity). This is only true in the Weyl basis that we are working in and will 
not be noticed in a non-Weyl basis. Less trivially, the correction term also vanishes for the de 


Sitter expansion (10), since now the two integrals in (35) compensate each other. Thus de Sitter 


is also an exact solution of the effective Friedmann EOM (35), although, again, this will not be 


true in a different (i.e. non-Weyl) basis for the non-local terms. 


Matter or dust instead is not an exact solution as the correction terms in (35) do not vanish 


for this case, and so we treat it iteratively. Since we are only interested in the LO corrections 
in this paper, we can treat them as perturbations above the unperturbed background and use 
the classical solution as an input to calculate them in an iterative manner. Following this, if we 


insert the matter Friedmann solution (12) in (35), we find 


H\t) - 


2567rGa 

3t4 


log/it + log 1^-1 


+ R1-1 

3 V«o 


SvrG 


(36) 


This iterative procedure converts the integro-differential equation into a simple differential equa¬ 
tion [IH]- We have assumed t > to, with to being a reference time, such that a{to) = 1 to solve 
this equation. Notice that, due to the weaker time dependence of the correction term, it remains 
subdominant with respect to the classical solution and therefore, will not affect the asymptotic 
evolution of the universe. The coefficient a, from Table is given by 


1 


a = 


(dvr) 




(37) 


where No is the number of scalars and x = 0(1) foi" rninimally (conformally) coupled fields. 
Note that for a reasonable number of scalars the graviton contribution to a is dominant. To 
obtain the numerical solutions of (36), we have set SttG = 1 and // = 1. In Figure]^ we have 
plotted the time evolution of the scale factor and the Hubble parameter for the classical solution 
and the LO corrections for the case of Nq = 100 minimally coupled scalars. As is evident from 
the figure, the magnitude of the LO corrections remains exceedingly small. With even larger 
number of scalars, one can expect that the magnitude of the LO corrections will increase but the 
character of the effect will not change. Moreover, as the scale factor evolves, the LO corrections 
die away due to their subdominant time dependence and the evolution is given by the usual 
classical solution. Although we have not attempted a complete numerical solution for this case, 
the analysis presented in the later sections indicate that the iterative solution indeed captures 
the subdominant effects. 


4.2 R^R gravity 

This correction term has been studied in the cosmological background, for the first time, in 
In this case, the non-local action is 


lNL = -j d^Xy/^R^R. 


(38) 


It is evident that this correction term represents a long range interaction between the scalar 
curvatures and is generally induced by all types of matter content, viz. scalars, fermions, photons 
and gravitons as clearly seen from TableTo derive the covariant EOM, we start with a variation 
of (38) with respect to the metric 
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5/nl = - 5^f^R-R + 2^/^6R-R - R-6n -R 
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'' (39) 

where we have used the exact relation = formally obtained by varying the identity 


□ i = 1. Using the variations (14) and (27) we immediately find the following correction to the 
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Figure 3: The time evolution of the scale factor (left) and the Hubble parameter (right) for the 
case of the logarithmic correction term in an expanding matter dominated universe with 100 
minimally coupled scalars (x = 0 ) and the graviton. 


Einstein’s tensor 
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167rG/3m2 


= 2G 
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□ 






(40) 


As a check, note that this relation, if evaluated in d = 2 (i.e. when Gfj_u = 0), gives exactly the 
(T^i/) of Polyakov’s action [53]. It also agrees with the correction term discussed in |22|. As done 


in the previous case, we can linearize the correction terms in (40) by introducing an auxiliary 
held, U = to hnd 


1 

IGvrG/Sm^ 


AG;.. 


-UU 


—‘2jGfii/U + 2gni,R + 2V;iV.t/ + ^ uU — aU 

R. 


It is now straightforward to write down the effective Friedmann EOM 


+ IdirG^m^ - 2HU - 2H^U^ 

U + 3HU 


SttG 

6(2H‘^ + H^ . 


(41) 

(42) 


Contrary to the case of the logarithmic correction term, the EOM in this case are a simple system 
of ordinary differential equations. Together with the continuity equation ([^, they form a closed 
system which can be solved for a given matter content of the theory. 

Radiation with ii = 0 is an exact solution of these EOM while de Sitter and matter are not. 
In these latter cases, we shall solve the equations (41) and (42) iteratively, as done earlier. For the 
de Sitter expansion given in (10), we hrst solve (42) to get the solution for U with the boundary 
conditions U{to) = U{to) = 0 , and insert it in (41) to obtain the correction term. The effective 
Friedmann EOM is then given by 


(^1 + {t - to) - e- 2 VM{t-to)^ ^ ^ ^ ( 43 ^ 

This iterative procedure reduces two equations to one effective equation which can now be easily 
solved for given initial conditions. The second term in the above equation arises due to the non¬ 
local correction which contains a constant term, a linearly growing term as well as an exponentially 
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Figure 4: The time evolution of the scale factor (left) and the Hubble parameter (right) for the 
case of correction term in a de Sitter universe with a single minimally coupled scalar 

(x = 0) and the graviton. 




Figure 5: Same plot as above but in an expanding matter dominated universe. 


decaying term in time. Note that, the correction term vanishes as t ^ to as well as A —>■ 0, as 
expected while at late times, i.e. t S> to, the exponentially decaying term can be safely neglected 
compared to the other two terms. It is also clear that the characteristic time scale involved is 
In order to solve this equation, we first need to read off the expression for (3 from Table 
[T] which is 


/3 = 


1 


(dvr) 


14 - 18x + 3x^ 
216 


13 


IVq T —iVi /2 —- -\~ — 

27 ins ^ 97 


108 


27 


(44) 


To obtain the numerical solution of (43) we have set SttG = 1, = 1 and A = 1 and we focus 

on scalars and graviton contributions to (3, although it can be easily done including fermions 
and photons. In Figure we have plotted the time evolution of the scale factor and the Hubble 
parameter for the classical solution and the LO corrections for the case of a single minimally 
coupled scalar (x = 0) and the graviton. It is evident from the figure that the LO corrections for 
the de Sitter case are significant. For the case of matter as an input solution we instead find the 
following effective Friedmann EOM 




1287rG f io\^ 
^7 ~ [j ) 


4 logy 



(45) 


which we have solved in the same manner and plotted the solutions for a and H in Figure As 
can be seen in this figure, the LO corrections for this case are slightly larger in magnitude than 
in the previous logarithmic case which can be understood due to the stronger time dependence 
of the non-local term in this case. 
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Finally, we would like to remark that it is possible to integrate the system (41) and (42) 


numerically to obtain the evolution of the scale factor and the Hubble parameter without the 
approximation involved in the iterative solution just exposed. But the outcome clearly shows 
that the iterative solution is indeed, not only more insightful analytically, but also quantitatively 
rather precise as we will show in section when discussing the equation of state parameter. 

4.3 R^R gravity 

This non-local correction is one of the most relevant IR term which has been argued as a simple 
and consistent non-local extension of GR and thus has been extensively studied in constructing 
models of dynamical dark energy [23 ED- For this case the non-local action is given by 


(46) 


^NL = j d'^x^/^R-^R. 


Note that, the effective action Inl is suppressed by a factor of 1/M^ compared to the Einstein- 
Hilbert term, as discussed in eq. (13) of [5]. From a dimensional point of view, it is easy to see 
that ^R^R behaves like a cosmological constant term and therefore, one might expect that 
it will play an important role in late time cosmology, as we will indeed see later. We also note 
that scalars and fermions, but not photons, generate this term. In order to obtain the covariant 


EOM, we start with a variation of (46) which leads to 


(5Inl = j d^x!^5V^R^R + 2^6R^R-2^R 




R 


(47) 


where we have used 5^ = — in the last term. By making use of the basic variations in 

(14) and (27), we obtain the following correction to Einstein’s tensor 




1 

□2^ 




1 


+ 9^iv^ ( i-i [ n2^ ] 


1 




(48) 


To express this equation in a simplified form, as done earlier, we introduce two auxiliary fields, 
U = r^R and S = ^R, and find 


167rG6m^ 


~ ‘2GnvS + 2guuU + 2VnV i/S' 


+ 2V^UV,S - 


-au = R 
-□s = u. 

One can now immediately write down these equations in the FRW background 

SvrG 


7/2 - 327rG<5m^ ( 2H'^S + HS + ^HS - ^US 

2 6 


U + 3HU 
S + 3HS 


-P 

Q[2H'^ + H 
U 


(49) 

(50) 

(51) 


(52) 

(53) 

(54) 


These coupled differential equations form a complete set which can be solved for a given back¬ 
ground evolution. It is now easy to rewrite these equations in terms of ' = ~ 

compare them directly with |31j . 
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Figure 6; The time evolution of the scale factor (left) and the Hubble parameter (right) for the 
case of R^R correction in de Sitter universe with a single minimally coupled scalar (y = 0) and 
the graviton. 




Figure 7: Same plot as above but in an expanding matter dominated universe. As shown in 
the inset of the figure on right, while the Hubble parameter for the classical solution behaves as 
the one induced by the LO corrections becomes nearly constant at late times, leading to an 
accelerated expansion of the universe. 


As in the earlier case, radiation with ii = 0 is an exact solution for this case while de Sitter 


and matter are not. We shall solve the closed set of equations (52), (53) and (54) iteratively for a 
given background expansion of the universe. As mentioned in the previous case, we first solve 


and (54) for the two scalar potentials U and S with the boundary conditions t/(to) = U{tQ) = 0, 
S{to) = S{to) = 0 and insert these solutions in (52) to obtain the explicit form of the correction 
term. On a de Sitter background, the effective Friedmann EOM reads 


- 


USttG 


27 


6 m^< 4 — 3\/^ (t — to) + 3A(t — to)^ — [s — (t — to) 


g-Vav (t-to) 


+ 


4 + 2\/M(t-to) 




A 


(55) 


The correction term in this case contains a constant term, a linear and quadratic term as well as 
exponentially decaying terms of characteristic timescale As earlier, the correction term 

vanishes as t ^ to as well as A —>■ 0, as expected. At late times, the exponentially decaying terms 
can be ignored compared to the others and thus, the most relevant term in the deep IR regime 
will be the quadratic one. As said, this non-local term is not generated by photons but only by 
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scalars and fermions and therefore, the expression for 6 is given by 


5 = 


1 

(47r)2 


15 — 6X — AT 1 AT 

^-iVo + -A^ 1/2 + 


144 


17 


(56) 


We now solve this equation numerically and plot the time evolution of the scale factor and 
the Hubble parameter for the classical solution and the LO corrections for the case of a single 
minimally coupled scalar (x = 0) and the graviton in Figure]^ We find that due to the stronger 
time dependence of the correction term, the LO corrections are significantly large at late times. 
To understand the effect of the correction term in the case of matter, we first write down the 
effective Friedmann EOM 




327rG , 4 f 119 

81 I 3 


26 log-- 

^ t t 





(57) 


and then solve it numerically. Note that, the correction term correctly vanishes as t —)• to- In 
Figure we have plotted the numerical solution for the scale factor and the Hubble parameter 
for the case of a single minimally coupled scalar (x = 0) and the graviton. We hnd that the LO 
corrections are quite significant for this case as well. More interestingly, we notice that, at late 
times, the Hubble parameter induced by the LO corrections becomes nearly constant and leads 
to an accelerated expansion of the universe, as shown clearly in the inset of the right Figure 

We shall further discuss this interesting behavior in the next section when discussing the 
equation of state parameter. We will also show that the iterative solution just exposed is, as in 
the previous case, quite a good approximation in comparison to the full numerical solution in 
both the cases. 

Before we close this section, an interesting remark about the corrections in de Sitter for ^ 
and terms is in order. It can be seen from (43) and (55) that the corrections due to the non¬ 
local terms in de Sitter lead to an effective time dependent cosmological constant in a consistent 
way which is a distinct outcome of our formalism and it is induced in a completely covariant set 
up without breaking the diffeomorphism invarianc^ This leads to interesting phenomenological 
consequences, as discussed in [5411^(5^l57] . 


5 A minimal comparison of non-local corrections 

In this section we quickly compare the iterative solutions obtained in section for the non-local 
corrections predicted by the EFT of gravity in order to understand their relative role in the 
evolution of the universe at late times. We can treat these correction terms as an effective dark 
energy component and write the general Friedmann EOM in the following simple form 

H^ = ^{p + Pi,). (58) 

The effective dark energy densitiy for each of the various LO non-local corrections can be 

^Note that the canonical method of calculating one-loop quantum corrections also does not break the diffeomor¬ 
phism symmetry. However, the canonical method is not manifestly diffeomorphism invariant while our covariant 
approach is. 
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Figure 8: The equation of state parameter t^de is plotted as a function of time corresponding to 
the different non-local correction terms. It is evident from the figure that the corrections induced 
by matter for the case of the R^R term can indeed lead to an accelerated expansion of the 
universe with tCde < — 1- Note that the behavior of rcde for different cases is independent of the 
respective coefficients a, /3 and 6. 


read off from the previous section and have the following form 
32a 


a,m /,\ 

Pde (*) 


( 5 , A /, \ 

Pde (*) 


5,m /,\ 

Pde (*) = 


16 


logmt + log ( ^ - ij + ^ - 1 


(l + {t - to) - e-2v^(*-*o)^ 


9 ^ 


l-(^) -41og^ 


(59) 

(60) 

(61) 


1 fi ( 

y4 — 3\/^ {t — to) + 3A(f — fo)^ 

- 8 - \/^ {t - to)l + [4 + 2\/^ {t - to) 




(62) 


Ifed -™-261og^2 + ^(y?- 

27 1 3 ^ t t \ 3 


8Iog‘f 


-8,‘f 


2 /_ 4 

3 I f / 3 I f 


(63) 


Here, the superscripts on each term indicate the corresponding correction term and the input 
classical solution. For instance, Pdg™ indicates the correction term for the logarithmic case with 
matter solution as an input. Note that, as also discussed earlier, the correction term for the 


logarithmic case vanishes for the de Sitter expansion (10) as all the radiation cases. The aim is 


to understand the late time evolution of these corrections and see if they can mimic an effective 
dark energy component, leading to an accelerated expansion of the universe. The fact that the 
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Figure 9: We compare the exact and iterative solutions for w^e for the R^R (left) and R^R 
(right) correction terms. We found that the iterative solution is a very good approximation to 
the exact solution in both the cases asymptotically. 


correction terms are covariantly conserved leads to the continuity equation for pde; 

Pde “b 3Ff(l -|- U?(Je)/^de — 0 i 


(64) 


where tCde is the effective equation of state parameter for the dark energy component. This 
equation immediately leads to 

- 3^ ()i) ■ 

It is interesting to note that while pde for different cases depends on the respective coefficients 
a, f3 and 5, the equation of state parameter Wde is indeed independent of them. One should 
typically plot the evolution of the total equation of state parameter w for a system with mnltiple 
fields or fluids but it turns out that, using the continuity equations for the total energy density, 
matter and the induced dark energy component, one can write down an expression for w as 


w = - —Wde , (66) 

Pra + Pde 

which tells us that w ^ Wde for pm ^ Pde at late times and therefore, we plot Wde as a function 
of time for the different cases, as shown in Figure While the corrections induced by matter 
in the cases of logarithmic and ^ terms lead to Wde > 0 and —1/3 < Wde < 0, respectively, it 
is a naive expectation that the correction terms induced by a de Sitter solution in any case will 
lead to Wde < which is what we have found for the cases of q and terms. Moreover, 
the corrections induced by matter in the case of ^ are indeed the most interesting ones as they 
generically lead to a phantom like behavior with Wde < — 1 and can indeed explain the current 
acceleration of the universe without the need for a quintessence field [281131] . 

Before we close this section, a few comments are in order about the comparison between the 
exact and iterative solutions for the cases of R^R and R^R corrections. We have plotted both 
the solutions for these two cases in Figure It is clear from the figure that the iterative solution 
is indeed a good approximation to the exact solution in both these cases and seem to capture the 
subdominant effects quite well. Therefore, we can argne that the iterative solution serves as a 
good approximation to the exact solution in all the cases. Hence, it is a very useful tool to obtain 
the solutions in cases wherein the system of equations can not be solved exactly. 


6 Discussion and Outlook 

In our previous paper [^, we had constructed an EFT of gravity in a complete covariant manner 
by means of an expansion in the inverse powers of Planck mass and identified the classical theory 
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together with the leading and next-to-leading quantum corrections. We then computed the 
leading order corrections to the quadratic order in curvature and found that these corrections 
in the effective action can be classified into two categories: local and non-local. In the FRW 
spacetime, the leading local correction is given by the E? term and its coupling is a free parameter 
which must be fixed by observations. The non-local sector instead contains different terms 
relevant at low energy and their coefficients are completely determined by the EFT of gravity 
once the matter content of the theory is specified. Since the E? term is a UV term and is primarily 
relevant during the early evolution of the universe, an inflationary scenario based on it can be 
constructed [SOj. As we discussed, this scenario can be naturally accommodated in framework of 
the EFT of gravity. 

In this paper, we have focussed our attention on studying the cosmological consequences of 
the non-local terms when each of them is combined with the Einstein-Hilbert action. We have 
explained how to derive the covariant EOM for all the cases and also how to correctly project 
these expressions on the ERW background, taking care in explaining how to represent the different 
non-local contributions. It turns out that these terms are not so difficult to handle and can be 
easily dealt with in our covariant formalism. We have then studied both the analytical and 
iterative solutions in all the cases and also commented on the comparison between the iterative 
and numerical solutions whenever possible. We have compared the evolution of the scale factor 
and the Hubble parameter induced by the non-local corrections to the classical solutions in all 
the cases and found that the imprints due to the E^E term have the strongest effect among all 
the terms considered. Einally, we compared all these contributions after rewriting them in terms 
of an effective dark energy component characterized by an equation of state parameter and found 
that E^E term can indeed drive an accelerated expansion of the universe at the present epoch 
[281 El]. 

One of the key result of this paper is that all the non-local terms in the EET of gravity arise in 
a consistent way from the first principles and have precise coefficients which depend ultimately on 
the matter content. Over the last few years, a large class of phenomenological models containing 
one or more of these non-local terms have been proposed in the literature. Here we have shown 
that all such terms appear with precise coefficients which hints at the possibility to falsify these 
phenomenological scenarios. Also, all these non-local terms generally appear together in the EFT 
of gravity unless certain matter content of the theory renders one of them non-vanishing while 
killing the others. But, in any case, the logarithmic term is always present and therefore, none 
of the other non-local terms can be considered alone. 

One of the opportunities offered by the covariant EFT of gravity is the inclusion of the 
cosmological constant in a consistent way. The combined effect of the cosmological constant and 
the non-local terms is to effectively induce a time-dependence of the former [SU ES ES] EZj- 
We emphasize that this is an interesting outcome in our covariant formalism which is manifestly 
diffeomorphism invariant contrary to the canonical method. We plan to come back to it in the 
future. 

An interesting possibility would be to study a scenario wherein the leading local correction 
term is combined with the IR relevant non-local terms. Such a scenario would lead to an early 
inflationary epoch together with an accelerating phase at the present epoch thereby unifying 
inflation and dark energy in a single framework by means of pure gravitational corrections without 
the usual need of a scalar field [58l ESI EQ] . Furthermore, the next natural step would be to study 
the evolution of cosmological perturbations in the EFT of gravity. As mentioned earlier, the 
Weyl contributions to the effective action will play a pivotal role in this context and it will be 
interesting to look for such imprints in cosmological observables. 
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A Non-local operators on FRW 


We are interested in providing a precise meaning to the Green’s function and the operator 

A.l Green’s function 

To compute the Green’s function on the FRW background, we need to solve 

(—□ + — t'] m) = 5{t — t') 

where 


□ = - 


1 d 
a^{t) dt 




and 


5{t -1') = 


6{t — t') 


(67) 

( 68 ) 

(69) 


is the covariant delta function, normalized such that f dt'a^{t')5{t — t') = 1. Plugging (68) and 
(1^ in ((67|) leads to 


1 d 
a?{t) dt 


K^’dt ' ’ 7 ^a?{t)a^{t') 


(70) 


The full kernel in a general FRW spacetime is not realistically tractable. So as an approximation 
we drop the time derivatives of the scale factor and rewrite (70) as 


[7(t)G^^^(t - t'; m)] + 0(a) = 5{t - t'). 


(71) 


Thus, for a slowly varying scale factor, we can use the retarded (or casual) flat space massive 
Green’s function 


07*(t;m) = 0(t) 


sin mt 


m 


to write the following symmetric form for the solution in €3 

a™{t - t';m) = ^(*-0 

■ 7a3(t)a3(F) rn 


(72) 


(73) 


In the massless limit we find 


Gl^^{t-t';0) 


e{t-t'){t-t') 

^/a^{t)a?{t') 


+ 0(d). 


The general FRW Green’s function can be found using the Strum-Liouville theory. 

A.2 The log ^ distribution 

The distribution = log 7^ is precisely defined as 


^FRW 



ds 


1 

+ s 


G^^^(t-t'-,^/^) 


(74) 


(75) 
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where different boundary conditions on the Green’s function different kernels . 

We are interested in the retarded kernel obtained by plugging (73) in We need to 


evaluate the following integral of the retarded flat space Green’s function 

roo i-oo • /-j. roo 


I(t)= / dsG_(t;y/s) = 6{t) / ds- _ 

Jo Jo 


= 29{t) / dy sinyt, 

Jo 


where we have changed variables to y = y/s, dy = The last integral is meaningful only if 
interpreted as a distribution and in particular we hav^ 


fOO ^ 

/ dysmyt = V-, 
lo ^ 


where the principal value distribution is defined as 




e ^-0 


" g{x) 


dx 


fix) 


Xi\f{xi)=0 °° 

Thus, bringing the step function inside the principal value, we get 


r 

x.+e fix) 


lit-t') = 2V- 


t -1' 


(76) 


(77) 


(78) 


Integrating against a test function gives 
dt' I{t-t')f{t') = 2V I 


f 




' —OO 


t - t' 


= 2 lim 

e-s>0 


= 2 lim 

e-S>0 




dt 


t - t' 
,9{t — t' — e) 


=0 


or simply 


V 


t-t' 

- t') _ -t' -e) 


fit '), 


t - t' 


= lim 


e^-0 t — t' 


The other integral in (75) needs a regularization 


I 


OO 

ds—F. -= lim / 

y + s e-^-O Jq 


ds—F, -= —2 lim log u e . 

yj + s e^O 


(79) 


(80) 


As a distribution this is represented as —25{t — t')logye and combining with (79) finally gives 
the flat space kernel 

9{t — t' — e) 


LyUt-t') = -2 lim 
e-S>0 


t - t' 


+ 5{t - t') log/re 


Using (|73|) we can write the following expression for the FRW retarded kernel 

9{t — t' — e) 1 


Ll^^{t-t') = -2 lim 
e-S-O 


5{t-f) 

H- , „, . „ , , = log y e 


■\/aJ{t)a^{t') t — t' yJa^{t)aJ{t 
which is valid for a slowly varying scale factor. 


(81) 


(82) 


^The calculation steps are as follows 


f 

Jo 


lim [ 


dy sin yt = lim j dye^^ sin yt=J lim 


1 


t + ie t — ie 


= r- 


where in the last step we used lime_»o ~ Tf T 
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